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Perturbation Theory on the Transition Temperature and Electronic 
Properties of Organic Superconductor 
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We study the superconducting transition temperature and the electronic properties of the 
metallic phase of K-type (BEDT-TTF)2X which shows unconventional properties in experiments, 
on the basis of the third order perturbation theory for a simple effective Hubbard model of a 
nearly triangular lattice. Appropriate transition temperatures and d^2_y2 symmetry of the 
gap function are obtained in good agreement with experimental results. We also calculate the 
transition temperature by the fluctuation-exchange approximation(FLEX) in order to compare 
the two approaches; FLEX gives higher transition temperatures rather than the perturbation 
approach. However, it is also found that the vertex corrections, which are ignored in FLEX, have 
a crucial effect on Tc for strongly frustrated systems. The density of states and the normal self- 
energy calculated in this perturbation scheme show the nature of the conventional Fermi liquid 
near the Mott-insulator. Thus, our perturbation approach is applicable to the conventional 
metallic phase of this compound, while it cannot explain the (pseudo-)spin gap phenomenon 
which signals the non-Fermi liquid. 

KEYWORDS: perturbation approach, electron correlation, spin fluctuation, vertex correction, Fermi liquid, d-wave 
superconductor, density of states, self-energy 



§1. Introduction 

Today, there are many organic compounds which 
show the superconductivity at low temperatures. One 
of the most interesting superconductors among these 
compounds is K-type (BEDT-TTF)2X (X=Cu(NCS)2, 
Cu[N(CN)2]Br, Cu[N(CN)2]Cl and k-(ET)2X for an ab- 
breviation). These compounds show superconductivity 
under the pressure appropriate to each compound. We 
have two ordered phases by varying pressure, one of 
which is antiferromagnetic insulator in the lower pressure 
region and the other is superconductivity in the higher 
pressure region. El 

By the nuclear magnetic resonance(NMR) experi- 
ments, it has been confirmed that, below the transi- 
tion temperature(Tc), there are no coherence peaks, the 
Knight shifts decrease and the spin lattice relaxation rate 
behaves as These results imply that the symmetry 

of the Cooper pair is not the isotropic s-wave, but the 
anisotropic symmetry and spin singlet, indicating the d- 
wave. These properties are quite similar to the cuprate 
superconductors (high-Tc)- The mechanism of this su- 
perconductivity cannot be explained on the basis of the 
conventional phonon-mediated superconductivity, which 
usually shows s-wave property and the isotope effect. 

Theoretically the unconventional superconductivity 
have been discussed keeping high-Tc cuprates in mind, 
roughly speaking from two approaches which are the thj 
ory based on the resonating- valence-bond (RV|B)|-statel 
{t-J model) and the spin-fluctuation model.lj^'B''la' In 
these approaches RVB theory assumes a strong on-site 
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Coulomb interaction, and is expected to be effective to 
the doped Mott insulator like the high-Tc materials in 
the underdoped region. On the other hand, our system is 
metallic although it has the nature of the half-filling; this 
fact suggests that the on-site Coulomb interaction is not 
so strong. Therefore, we consider the weak-coupling ap- 
proach is appropriate to this superconductor. When we 
discuss the spin-fluctuation mediated superconductivity, 
we have several methods of determining the attractive 
force. Moriya, et al. determined Imx('7,w) (the imag- 
inary part of the dynamical susceptibility) by the self- 
consistent renormalization(SCR)p' and Pipes, et al. de- 
termined it by NMR experimental results. B They solved 
the Eliashberg equation by using the phenomenologically 
determined susceptibility, and obtained the appropriate 
Tc- Before them, simple calculations by the Random 
Phase Approximation (RPA) were done by Shimahara et 
aln However, the Tc obtained by this method was rather 
lower than those of the above methods owing to thejOver- 
estimated damping effect as pointed out by Hotta.EP (It 
should be noted that Shimahara also explored the possi- 
bility of spin-fluctuation mediated superconductivity in 
the organic superconductor (TMTSF)2X.) 

These theories are expected to be effective to k- 
(ET)2X which shows strong antiferromagnetic spin fluc- 
tuation as observed by NMR experiments. Here, in order 
to understand these compounds on the basis of a more 
microscopic model, we take the following approach. 

Because this superconductivity occurs only near the 
Mott-insulating phase under the controlled pressure, the 
electron correlation on the same site is expected to 
play an important role on the appearance of this su- 
perconductivity. One of methods to investigate elec- 
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tron correlation-induced superconductivity is the per- 
turbation method, with respect to Coulomb repulsion 
U. In the studies of the transition temperature for the 
high-Tc cuprates, this method was adopted by HottaE^ 
He calculated Tc by the third order perturbation treat- 
ment(TOPT), and discussed the effect of the normal self- 
energy and the vertex correction on the calculated 
and then found out that the vertex corrections suppress 
Tc, compared with that obtained by a summation of the 
RPA-like diagrams. 

Another method of investigating this kind of supercon- 
ductivity from the weak coupling regime is the fluctua- 
tion exchange approximation(FLEX)j3'' which is a kind 
of the self-consistent calculation based on the RPA-like 
diagrams. The use of this method is reasonable because 
the anisotropic superconductivities usually appear near 
the magnetic ordered phase by controlling the pressure 
or the doping of carriers. Therefore, the summation of 
the specific diagrams which are favorable to the spin fluc- 
tuation may be justified. However, in this work we con- 
centrate on TOPT, without assuming the strong spin 
fluctuation like FLEX because it is not clear yet whether 
taking only the effect of particular diagrams is reliable 
for revealing the mechanism of superconductivity. 

The perturbation approach is sensitive to the disper- 
sion of the bare energy band by its nature, it implies 
that the lattice structures and the band filling play the 
essential roles in the calculation of Tc and the electronic 
properties. Although the TOPT calculation is already 
performed for the high-Tc, it is important for the study 
of this compound to calculate Tc newly on the basis of 
the effective model of k-(ET)2X. 

In this paper we calculate Tc for the simple effective 
model by TOPT and investigate its ?7-dependence and 
the dependence on the effect of the frustration. Re- 
cently the calculations of Tc by FLEX were performed 
by Schmaliaro^ and Kino et at, Kondo et ali3 on the 
two-band model and the single-band model respectively. 
Therefore we also calculate Tc by FLEX because these 
calculations were performed separately and there is no 
comparison with other methods yet; it is meaningful 
to compare the results by TOPT and FLEX. We also 
estimate the contribution of the vertex corrections, by 
comparing the results of the full TOPT calculation with 
that of only RPA-like terms included in anomalous self- 
energies. 

Another prominent property of this compound is the 
appearance of the (pseudo-)spin gap, which was_revealed 
by the NMR experiments in the metallic phaseO There- 
fore, it is necessary to study the electronic properties in 
order to see what extent our model and approximation 
which are used to calculate Tc are pertinent to. We cal- 
culate the density of states and the normal self-energy, as 
the physical quantities characterizing the metallic phase 
of this compound, and then discuss the results and fur- 
ther problems. 

§2. Formulation 

2. 1 Hamiltonian 

First we define the effective Hamiltonian in a minimum 
form which describes k-(ET)2X for the calculation of Tc. 



The Fermi surface(FS) of k-(ET)2X was studied by the 
Shubnikov-de Haas effect and found that it is well repro- 
duced by the tight binding appr|aximation based on the 
extended Hiickel approximation^' However it is known 
that a pair of ET molecules is considered as the basic 
structural unit (which is called a dimer) of the conduction 
sheet of k-(ET)2X because the intradimex-hpoping term 
is more than twice of the interdimer oncy'liJ Here, we 
briefly discuss the dimer model, which can be described 
as Fig. 1. 
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Fig. 1. (a) The original lattice structure of ft-(ET)2X when 
BEDT-TTF molecules are written explicitly. The line represents 
each molecule, f^j is the intradimer transfer integral, ti , t2 and 
t3 constitute the interdimer transfer integral, (b) The lattice 
structure of k-(ET)2X when the pair of molecules is dimerized. 
The circle represents a dimer, and the area which are surrounded 
by dashed lines, represent the unit cell in this case, t and t' are 
transfer integrals connecting the sites. 



The arrangement of BEDT-TTF molecules are shown 
in Fig. 1 (a). When a dimer is considered as a structural 
unit, the energy splitting between the bonding and anti- 
bonding orbitals is approximately given by 2i^. There- 
fore we can neglect the mixing between the bonding 
and antibonding orbitals, and the dimer model is ver- 
ified. Focusing on the antibonding orbital, each dimer 
is connected to the nearest sites by two kinds of trans- 
fer integrals. These transfer integrals can be written as 
t = {ti + ^2)72 and t' = t-ijl. (The transfer integrals ti, 
ti and ^3 correspond to those of Fig. 1 (a).) Thus, the 
lattice structure can be described as Fig. 1 (b) by using 
the above value of t and t' . 

Strictly speaking, transfer integrals t in Fig. 1 have 
two different values which cause the gap between the 
quasi-lD FS and quasi-2D FS and put two dimers in the 
unit cell. However, since the diffemnce between them is 
quite small (only 2 % of its value), Q we consider its dif- 
ference can be neglected and regard the minimum model 
as the lattice with one site in the unit cell. 

For the reasons stated above we adopt the following 
Hubbard Hamiltonian 

<zj>,cr <i,k>,a 

-I-J7 ^ n^^ini^X -f- /i^Tii^o-, (2.1) 

i <7 

where <j is spin indices, < i,j > indicates taking summa- 
tion over nearest neighbor sites and < i,k > over next 
nearest sites only in one direction as in Fig. 1. /i is the 
chemical potential which is determined so as to fix the 
electron number to 1 per site. The noninteracting energy 
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dispersion is; 

— — 2t(cosfcj. + cosky) — 2t'cos{kx + ky). (2.2) 

According to the quantum chemistry calculations^' we 
take t' /t ~ 0.7 as the realistic value of k-(ET)2X. It 
should be noted that this lattice has strong frustration 
and can be called a nearly triangular lattice. 

2.2 Green's function and self-energy 

The non-interacting Green's function is the following; 

Go(fc,e„)==: i (2.3) 

where e„ = 7rT(2n + 1) (n is integer) is the fermion- 
Matsubara frequency. The diagrams of the normal self- 
energy are shown up to the third order of interaction 
in Fig. 2 (the Hartree term is included in the shift of 
chemical potential). 
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Fig. 3. The diagrams of the RPA parts{which are included in 
FLEX) of the anomalous self-energy up to the third order. 



EPS File fig4.ps not found 



Fig. 4. The diagrams of the vertex correction parts(which are not 
included in FLEX) of the anomalous self-energy up to the third 
order. 



Fig. 2. The diagrams of the normal self-energy up to the third 
order. The solid and dashed lines correspond to the bare Green's 
function and the interaction, respectively. 



The analytic form of the normal self-energy is given 



by 



Sn(fc, e«) = T7 X! [^^^o(fe - fc', en - £„') 



fc', 



+U^4>l{k + fc', e„ + e„')]Go(fe', £«')• (2-4) 
In this equation xo ^-nd (pa are represented as 



k.n 



(2.5) 



B(q,cj,„) = -— ^Go(fc, e„)G'o(g — fc,^m - £«)■ 



N 



fc. 



(2.6) 

The anomalous self-energy can be divided into two parts, 
one of which is included in FLEX, and the other is not 
included in FLEX, i.e. the vertex correction. These dia- 
grams are shown in Fig. 3 and Fig. 4, respectively. 



Analytically these diagrams are written as following 
equations. 

Sa(fc, e„) = Sj^pA(fc, e„) + SyertC^' ^n), (2.7) 



where A:=(k,e„). 

2.3 The particle number 

The particle number is 1.0 per site in the real sys- 
tem. However, in our perturbation scheme, there is 
a discrepancy between the particle number based on 
the bare Green's function (which is denoted by no) ^^nd 
that based on the dressed Green's function(which is de- 
noted by n) in our perturbation scheme. The bare 
susceptibility(xo(g, i^m)) plays an important role in the 
calculation of Tc, that is, it determines the magnitude 
and the spatial and temporal variation of the interac- 
tion between the particles. To incorporate the nature of 
half- filling, we put no = 1, and require the number con- 
servation between hq and n. To satisfy this requirement 
we introduce the shift of chemical potential Sfi, which is 
included in the following dressed Green's function; 

G(fc, e„) = -. \ (2.10) 

and the particle number is given by 



= 2^^G(fc,6„). 



N 



(2.11) 



To satisfy n — 1.0 up to the third order of the interaction, 
expanding eq. (2.11) with regard to — I]n(fc, £«) , and 
using no = ^Y^k.n Go{k, e„) = 1.0, then we get 



(5/i = — 



N 



Efc,„G'§(fc,e„)Sn(fc,e„ 



Xo(o,0) 



(2.12) 
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SRPA(fe, e«) = J2^U + C/\o(fc + k') + 2U\lik + k')]F{k'), (2.8) 
fc' 

Svert(fc> en) = ~U^^ J2 Go(fc')(xo(fc + k') ~ Mk + k'))Go{k + k' - k^)F{k^) 

k' ,ki 

-U^^ E Go(fc')(xo(-fc + fc') - M-k + fc'))Go(-fc + fc' - fci)F(fci), (2.9) 



2.4 Eliashberg equation 

For the calculation of Tc, the Dyson-Gor'kov equation 
is linearized as 

Fik, en) = |G(fe, e„)|'Sa(fc, £„)• (2.13) 

Then the Eliashberg equation is given by 



When the eigen value of this equation is 1, the system 
is considered to be superconducting i.e. T = Tc- We 
don't assume the symmetry of the Cooper pair in calcu- 
lating Tc, which differs from Hotta,!]'' but calculating Tc 
and determining the symmetry are performed simulta- 
neously. 

§3. Calculated Results 

3. 1 Details of the numerical calculation 

We take the transfer integral t as the unit of the en- 
ergy arui,put t = 1. From the quantum chemistry calcu- 
lations Jlj' f ~ 70 ~ 80 meV, then Tc which is experimen- 
tally measured around lOK, is normalized to Tc = 0.011. 
The full bandwidth is from 8 to 9 when the value of t' 
varies from to 1. 

To solve the Eliashberg equation we have to calculate 
the summations over the momentum and the frequency 
space. Since all these summations are in the convolution 
forms, we can carry out these by using the algorithm of 
the Fast Fourier Transformation. For the frequency, ir- 
respective of the temperature, we take 1024 Matsubara 
frequencies, which makes the frequency cut-off compa- 
rable to the bandwidth at T ~ 0.0020. On the other 
hand, for the momentum space, in cpier to avoid the fi- 
nite size effect at low temperaturesjia^ we set points as 
fine as possible and divide the first Brillouin zone into 
128 X 128 meshes. 

3.2 Transition temperature 

We calculate Tc by solving the Eliashberg equation 
(2.14). The gap function shows the node at kx = ky and 
kx = —ky and changes the sign across the node in all 
approximations and for all parameters, namely 

Sa(fc, En) oc cosfca; — cos/cy. (3-1) 

The symmetry of Cooper pair is dx^^y-i. The U- 
dependence of Tc is shown in Fig. 5. 

In this figure, we show also the results obtained by 
FLEX calculation and that for only the RPA-like dia- 
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Fig. 5. The calculated Tc. TOPT(RPA-only) in this figure means 
that only RPA-like diagrams of anomalous self-energies up to 
third order arc included. (All normal self-energies are included 
up to the third order.) The diagonal transfer t' = 0.7. The unit 
of energy is the transfer t. 



grams(of anomalous self-energies) included in TOFT, in 
addition to TOFT calculation for comparison. This re- 
sults indicate that for larger U higher Tc are obtained 
commonly for all approximations of calculations. This 
is quite natural for our perturbation scheme and corre- 
sponds well to the experimental-results which are mea- 
sured by applying the pressure;E-J the pressure increases 
the bandwidth W and then U /W is reduced. 

We can see the differences between TOPT and FLEX 
in Fig. 5; Tc calculated by FLEX is higher than TOPT 
for moderate values of U. This is because RPA-like di- 
agrams are included up to higher order in FLEX, and 
the spin-fluctuation is largely enhanced. This reasoning 
is supported by comparison between TOPT and TOPT 
(RPA-only); the RPA-like diagrams of anomalous self- 
energies(eq. (2.8)) are responsible for attractive interac- 
tion, while the vertex corrections (eq. (2.9), this term is 
not included in FLEX as noted above) are responsible for 
repulsive interaction and then have the effect^of lowering 
Tc (This was firstly pointed out by Hotta.Q'). There- 
fore it can be said that the mechanism based on TOPT 
is the spin fiuctuation induced one. The another differ- 
ence between the results of TOPT and FLEX is their 
[/-dependences. The rate of increase of Tc as U is varied 
in FLEX is shghtly smaller than that of TOPT. This is 
probably due to the fact that the damping rate in FLEX 
is large because it is largely enhanced by the RPA dia- 
grams, as well as the attractive interaction is enhanced. 

There are maiui, quantum chemistry calculations as 
mentioned aboveEj^ and the obtained values for the 
transfer integral t' are roughly from 0.5 to 0.8, but it 
is difficult to know the real value of the transfer integral. 
Therefore it is meaningful to study the i'-dependence of 
Tc. The calculated results are shown in Fig. 6. 
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Sa(fc) = ^[C/ + U^xoik + k') + 2U\l{k + fc')]|G(fc')pSa(fc') 

k' 

-f^'^ Go(fc')(xo(fc + k') - Mk + fc'))Go(fc + fc' - fci)|G(fci)|2Sa(fci) 

E Go(fc')(xo(-fc + fc') - <^o(-fc + fc'))Go(-fc + fc' - A:i)|G(fci)pEa(fci). (2.14) 

k' 
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Fig. 6. The calculated t'-dependence of Tc in TOPT. The value 
of U is fixed to 6.50. 



The obtained highest value is Tc = 0.06400 for t' = 
0.30. For this value oiU — 6.5, and within our numerical 
accuracy, is not obtained for t' > 0.8. To obtain 
Tc for this region of t' , it is needed to increase U . At 
t' = 0.00, the Fermi surface is of perfect nesting and the 
bare susceptibility diverges. For this case we couldn't 
obtain the Tc. 

The calculated results of the static bare susceptibility 
(eq. (2.5)) are shown in Fig. 7 for various values of t' . 
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Fig. 7. The momentum dependence of the static bare suscepti- 
bility for various values of t' . These results are obtained for 
T = 0.0500. 



For the typical value t' = 0.70, the susceptibility has 
the incommensurate peak around (tt, tt) and asymmetric 
peak around (— 7r,7r). This result indicates the d^2_y2 
symmetry of the gap function. This momentum depen- 
dence of the haxe susceptibility is similar to the result 
of Kondo, et alLli> obtained by FLEX. Therefore it can 
be said that both of the approaches which are TOPT 
and FLEX use the nesting property to acquire the spin 
fluctuation and the attractive interaction. Here, the t'- 
dependence has the following feature. When t' is de- 
creased, the peaks at this momenta become sharper and 
then at t' = 0.10 the incommensurability and the asym- 
metry of the peaks disappear. These facts explain the 
results shown in Fig. 6, which generally shows higher 
Tc for smaller t' . Physically this corresponds to the fact 
that for smaller t' , the lattice is closer to the square lat- 



tice and the antiferromagnetic fluctuation is larger due 
to reduced frustration. Tc at t' = 0.1 is lower than that 
of t' = 0.3 in spite of its sharper peak for t' = 0.1. This is 
because as f decrease, the nesting property is enhanced 
and then the spectral density in the vicinity of the Fermi 
level moves to the two incoherent parts(i.e. which be- 
come the higher and the lower Hubbard bands in the 
Mott insulator). This transformation of the spectrum is 
shown in Fig. 8 
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Fig. 8. The density of states as t' is varied, at [/ = 6.50 and 
T = 0.06500. The inset shows this figure focused near the Fermi 
level. 



In order to see how the vertex corrections influence Tc 
when t' is varied, we also calculate Tc by including only 
RPA-like diagrams(Fig. 3) in anomalous self-energies, 
in other words, neglecting the vertex corrections. The 
results are shown in Fig. 9. 
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Fig. 9. t'-dependences of Tc, calculated by including only RPA- 
like terms in anomalous self-energies(All normal self-energies up 
to the third order are included.), and by TOPT. The value of U 
is fixed to 6.50. 



For a comparison with TOPT in the order of magni- 
tude, we take the linear-log plot as is shown in the figure. 
For the case that only RPA-like anomalous self-energies 
are included, we can also see that for smaller t' , higher 
Tc is obtained for the same reason as the calculation by 
TOPT(Fig. 6). However, as t' is increased, the rate of 
decrease of Tc by TOPT (which includes the vertex cor- 
rections) is larger than TOPT(RPA-only). This property 
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of t'-dependence is clearly shown in Fig. 10. 
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Fig. 10. i'-dependences of [Tc - T™="' [ /T™="' , these values are 
calculated based on Fig. 9. Here, T™'^ is the maximum value 
of Tc by each method(i.e. TOPT or TOPT(RPA-only)), when 
t' is varied. The value of U is fixed to 6.50. 



From Fig. 9 and Fig. 10, we can clearly see that 
for larger t' , Tc is largely reduced by taking the vertex 
corrections. At t' ~ 0.3, 0.4, the value of by TOPT 
is half or one third of that by TOPT(RPA-only). This 
ratio is comparable to the high-Tc calculation by HottaE^ 
However, at t' = 0.70 corresponding to the real system, 
the value of Tc by TOPT is reduced by one order of 
magnitude compared with TOPT(RPA-only). This fact 
suggests that for large t' , i.e. when the frustration is 
large like our system, neglecting the vertex corrections 
is not a good approximation. Although the reason is 
not so simple from the complicated form of eq.(2.9), this 
property is probably due to the fact that the peaks of 
bare susceptibility are not so sharp for larger t' as is 
shown in Fig. 7. Because Re0o(g, 0) doesn't have large 
t'-dependence as is shown in Fig. 11, xo{q, 0)— Re0o(9) 0) 
has more remarkable peak around (0, 0) for large t' as is 
shown in Fig. 12. 
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Fig. 11. The momentum dependence of Re</)o((jr, 0) from eq. (2.6) 
for various values of t' . These results are obtained for T = 
0.0500. 
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Fig. 12. The momentum dependence of xoC*?: 0) — Re</<()(q, 0) for 
various values of t' . These results are obtained for T = 0.0500. 



From Fig. 12, it is seen that for small t' , the peak 
around (tt, tt) makes the whole structure featureless. 



However, for large t', the strong frustration makes the 
peak around (0, 0) remarkable. Therefore the vertex 
corrections affect the gap function more repulsively for 
the d-wave symmetry at large t' from eq. (2.9). In 
other words, they suppress the antifcrromagnetic spin- 
fluctuation. Although within the third order perturba- 
tion, it can be said that the vertex corrections have a 
crucial effect on the calculation of Tc for strongly frus- 
trated systems. 

3.3 Density of states and self-energy 
The density of states(DOS) is given by 

p(uj) -- V ImG(fc, uj). (3.2) 

TT ^-^ 

k 

We calculate lmG{kjMj) from eq. (2.10) by using the 
Pade approximation.IlM) The f7-dependence of DOS is 
shown in Fig. 13. 
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Fig. 13. DOS at r = 0.01500, i' = 0.70 and various values of U 
as shown in the figure. The arrows indicate the incoherent peaks 
for each value of U. 



From this figure we can see that the incoherent peaks, 
which are corresponding to the upper and the lower Hub- 
bard peaks, grow as U increases. (The upper peak can be 
seen for small U, while the lower peak can be seen only 
for large U. This tendency reflects the original asym- 
metric structure of DOS.) This is the expected behavior 
for our perturbation scheme because for larger U we ap- 
proach the Mott-insulator phase, and coincides with the 
investigation on the infinite dimensions This behavior 
cannot be obtained within FLEX as is shown in Fig. 14, 
although U is very large and the system is expected to 
be close to the Mott insulator. 



EPS File figl4.ps not found 



Fig. 14. DOS calculated by FLEX at T = 0.02000, t' = 0.70 and 
various values of U as shown in the figure. 



The origin is probably ascribed to the self-consistent 
requirement in FLEX. In addition to this we see the 
asymmetric structure of DOS as expected for the nearly 
triangular lattice in Fig. 13 (i.e. because of the particle- 
hole asymmetry.). There is the coherence quasi-particle 
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peaks near the Fermi levels for all values of U, and we 
cannot obtain the pseudogap behavior within our model 
and approximation. 

The normal self-energy ia.obtained from eq. (2.4) by 
using Fade approximation.lia' The real part(ReI]„(fc, uj)) 
and the imaginary part(lml]„(fc, uj)) of the self-energy at 
the Fermi momentum are shown in Fig. 15 and Fig. 16 
respectively. 
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Fig. 15. The real part of the normal self-energy at the Fermi 
momentum ,at T = 0.01500, t' = 0.70 and various values of U 
as shown in the figure. 



EPS File figl6.ps not found 



Fig. 16. The imaginary part of the normal self-energy at the 
Fermi momentum ,at T = 0.01500, t' = 0.70 and various val- 
ues of U as shown in the figure 



The w-dependence of both parts near oj = are re- 
spectively given by ReS„(fc,Ci;) oc —uj and ImS„(fc,Ci;) oc 
— w^. (The small dip of the imaginary part at C/ = 7.0 
in Fig. 16, might be ascribed to the lack of accuracy of 
the Fade approximation.) This behavior is the expected 
one for the usual Fcmii liquid, and is opposite to the 
result of Maly, et a/.Ej-' which shows in their research in 



1 . , dRes„(k,u 



> u aiiu a i 

uj—0 — 

obtained by using f-matrix approximation, resulting m 
the appearance of the pseudogap behavior. Other points 
which should be noted in Fig. 15 and Fig. 16 are that, 
as U increases, the slope of ReI]„(feF,'^) 'At uj — be- 
comes steeper and the coefficient of cj^ for the imag- 
inary part becomes bigger. This indicates that the 
mass and the damping rate of the quasi-particle become 
larger (physically this corresponds to the increase of the 
resistivity) as U increases because the mass enhancement 
3ReS„(fc,i^) ^j^g damping 



and y ' 



> 



j=0 



0=0 



factor is given by 1 

rate is given by — lTOS„(fc, oj). These results are the typ- 
ical Fermi liquid ones, and cannot show the anomalous 
behavior like the pseudogap. 

§4. Summary and Discussion 

In this paper we have calculated Tc for the simple 
effective model of k-(ET)2X by solving the Eliashberg 



equation on the basis of TOPT, and have investigated 
the possible mechanism by the spin fluctuation on this 
superconductor. We have adopted the Hubbard model 
with a nearly triangler lattice in which one of these trans- 
fer integrals, t' — 0.7; we have taken J7 as a parameter 
corresponding to the experiment in which the pressure is 
applied, and fixed the system to half-filling in agreement 
with the real compound. The obtained results roughly 
agree with the experiment in the high pressure region in 
its [/-dependence and values. 

The results obtained here show that the superconduc- 
tivity in this organic conductor possibly arises from the 
same mechanism, by which we mean the antiferromag- 
netic spin fluctuation induced superconductivity, as that 
of the cuprate system, since both of them are connected 
by changing a single parameter t' . 

We have also compared the results obtained by TOPT 
and FLEX. For the moderate value U, of FLEX is 
higher than that of TOPT. This is because FLEX in- 
volves the higher order spin fluctuation in spite of its 
large damping rate. However, from our perturbative as- 
pect, because of the offset between large attractive inter- 
action and large damping which are both given by RPA- 
like diagrams, the appropriate values of Tc are obtained 
by FLEX. Therefore it should be noted that whether 
so large enhancement is physical one or not is a future 
problem. Concerning this point, we have shown in Fig. 
5 and Fig. 9 that the vertex corrections reduce Tc by 
one order of magnitude. Up to the third order terms, it 
is found that the vertex corrections reduce Tc of strongly 
frustrated systems like our system, more seriously than 
that of not so frustrated systems like high-Tc, although 
it cannot be said generally beyond the third order from 
our scheme. This fact suggests that the calculations of 
Tc which include only the RPA-like terms are question- 
able and should be carefully performed with the vertex 
corrections 

To investigate to what extent this scheme is effective, 
we study DOS and the normal self-energy. The obtained 
behaviors of both quantities are the expected ones for 
the Fermi liquid near the Mott-insulating phase. Espe- 
cially the upper and the lower Hubbard bands in DOS 
are obtained in our perturbation scheme, while they are 
not obtained by FLEX. From this behavior of DOS, it is 
doubtful that the physical quantities obtained by FLEX 
have the nature of being near the Mott-insulating phase. 
It can be said that one of the elements which improve 
FLEX is possibly to take this property into account. 

From the results of DOS and the self-energy, this per- 
turbation approach is valid in the conventional metallic 
phase in the phase diagram of McKenzieH^ However the 
(pseudo-)spin gap which is one of the most interesting 
properties in the strongly correlated electron system, is 
not reproduced in our scheme. 

From the analytic property in the power series expan- 
sion of U for the Fermi liquid, we cannot derive this 
unconventional property by using the present method 
like TOPT or FLEX. There must be the crossover which 
shows the breakdown of the Fermi liquid at a certain 
value of U as is already found in experiments. It has not 
been known yet that this breakdown is caused by the su- 
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perconducting fluctuation ^)i£3 or the antiferromagnetic 
long(or short)-range orderO^'cS or other mechanism. It is 
important to decide the critical term which characterizes 
the crossover and its effective region of the temperature 
and the magnitude of correlation. This is the problem to 
be solved to understand this compound and the nature 
of the non-Fermi liquid adjacent to the Fermi liquid. 
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